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We report measurements of the two-magnon states in a dimerized antiferromagnetic chain mate- 
rial, copper nitrate (Cu(N03)2 • 2.5D2O). Using inelastic neutron scattering we have measured the 
one and two magnon excitation spectra in a large single crystal. The data are in excellent agree- 
ment with a perturbative expansion of the alternating Heisenberg Hamiltonian from the strongly 
dimerized limit. The expansion predicts a two-magnon bound state for q ~ (2n -f l)nd which is 
consistent with the neutron scattering data. 



PACS numbers: 75.10.Jm, 75.40.Gb, 78.70.Nx 

INTRODUCTION 

The basic physics of the elementary one-magnon exci- 
tations of lower-dimensional quantum antiferromagnets 
can now be regarded as well established, both theo- 
retically and experimentally through studies of materi- 
als that are reasonably accurate realizations of the spin 
Hamiltonians. In contrast, higher excitations such as 
multimagnon continua and bound states have attracted 
relatively little attention. This topic may prove to be 
a fascinating area for the application of few- and many- 
body techniques, and will involve interesting and non- 
trivial results in band structure, band mixing, bound 
state formation, phase transitions through the formation 
of condensates of magnetic excitations ||, ||, |[ ||, ^, Q , 
and other collective phenomena. 

Although some aspects of the physics of low-lying mul- 
timagnon states can be inferred from model Hamiltonians 
using standard theoretical techniques, few experimental 
studies of these higher excitations have been reported to 
date. Experimental difficulties that have precluded such 
work include relatively weak couplings of probes to these 
higher excitations, dominant contributions from the low- 
est one-magnon excitations, and resolution requirements 
in energy and wavenumber that are beyond the capabil- 
ities of most techniques. 

High-resolution inelastic neutron scattering should 
prove to be an ideal technique for observing some of these 
higher magnetic excitations. One can control both en- 
ergy and momentum transfer, so that the existence and 
spectral weight of higher magnetic excitations can be es- 
tablished and quantified. The new generation of high- 



intensity neutron sources combined with high-resolution 
detectors should allow the observation of details of the 
multimagnon excitation spectrum such as band bound- 
aries, quantitative determination of the dynamical corre- 
lation function iS(Q, lo) and discontinuities within a band 
and weakly bound states just below the band edge. 
The principle limitation in this approach may be the un- 
avoidable A^* = 1 selection rule of magnetic neutron scat- 
tering, so that one can only reach spin S = 1 excitations 
given an S* = ground state and an isotropic spin Hamil- 
tonian. Other techniques such as Raman scattering can 
be used to study certain of these higher excitations, al- 
beit with strong constraints on the accessible spin and 
momentum quantum numbers. 

Systems that appear especially interesting for studies 
of higher magnetic excitations at present are quasi-lD 
spin chains and spin ladders, since many of these have 
gaps and hence will have separated bands of higher ex- 
citations and perhaps bound states of magnons. The 
alternating Heisenberg antiferromagnetic chain (AHC) 
with spin-1/2 is an example of such a system; with any 
amount of alternation < a < 1 (where a = J2/J1) 
the AHC has an energy gap in its one-magnon disper- 
sion with a second gap to the multimagnon continuum, 
and two-magnon bound states with spin-0 and spin-1 are 
predicted . The AHC is also attractive because of its 
relative simplicity and because of recent extensive theo- 
retical studies of the low-energy excitations in this model. 

In this paper we present results from an inelastic neu- 
tron scattering study of higher magnetic excitations in 
copper nitrate (Cu(N03)2 • 2.5D2O), which was recently 
confirmed by neutron scattering pC| ] to be an accurate 
realization of a strongly alternating Heisenberg antifer- 
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FIG. 2: Alternating chain layout showing dimers coupled to- 
gether. Each dimer is labeled by an index m, p is the sepa- 
ration between dimer spins, and d is the chain repeat vector. 

This model has a nontrivial spin-0 ground state, and 
for all a in the allowed range < a < 1 has a gap 
to the first excitation, which is a band of spin-1 excita- 
tions (magnons). In the "strong-coupling limit" a ^ 1 
the ground state approaches a system of uncoupled spin- 
dimers, and the one-magnon excitations can be accu- 
rately described as a single dimer excited to spin-1 (an 
"exciton"), delocalized on the chain to give states of def- 
inite along-chain wave number k. The energies and some 
matrix elements of these states have been evaluated as 
power series in the intradimer coupling a | pd| , p^ . 

To interpret experimental results approximate analytic 
forms of the wavefunctions and energies are useful and we 
include calculations of wavefunctions expanded around 
the single dimer eigenstates in the Appendix. The one- 
magnon S = 1 triplet has a gap energy A = J — a J/2 
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and dispersion to 0{a) 

'^imag(fc) = J — a J/2 cos(fc(i). (2) 

This magnon wavefunction can be visualized as a local- 
ized wavepacket of magnetic polarization (excitons) along 
x,y, or z carrying total spin-1 travelling through a fea- 
tureless singlet background with the gap energy coming 
from the effort expended in breaking a dimer bond. 

Approximate wavefunctions for the two-magnon states 
are given in the Appendix. At large separations the 
magnons do not overlap and they behave as free par- 
ticles, however when close they interfere and scatter off 
each other. In a one dimensional geometry such scatter- 
ing conserves particle number, energy, and momentum up 
to a lattice wavevector, and the scattering introduces a 
momentum dependent phase shift in the scattered wave- 
function. The energy of these states (to order 0(1/N)) 
is given by aJkiMW = wii„ag(fci) + Wimag(fc2), where 
fc = fci -|- fc2 is the total wavevector. The resulting con- 
tinuum using equation (2) is illustrated in Fig. 1 ^. 

As well as elastic scattering of magnons, bound states 
also form. To appreciate their physical origin consider 
two dimers, both in excited states, coupled by a sin- 
gle interdimer coupling a J. In the absence of coupling 
all the double excited states S — 0,1, and 2 have the 
same energy 2 J. However, the interdimer coupling splits 
these states. There is an 5 = 2 quintuplet of energy 
2 J -|- a J/4 (this energy is higher because the interdimer 
coupling favours antiferromagnetism whereas the S = 2 
states have all spins along the same direction - ferro- 
magnetic). There is also an 5 = 1 triplet of energy 
2 J — aJ/4, this lowering of energy is purely quantum 
mechanical and comes from resonance between the two 
excited dimers. Finally there is an S* = singlet of even 
lower energy, 2 J — a J/ 2, which gains resonance and also 
antiferromagnetic energy due to the spins in neighboring 
dimers pointing in opposite directions. 

As the S = and 1 states with excited dimers neigh- 
boring each other have lower energy than two well sep- 
arated excited dimers (by a J/4 and a J/2 respectively) 
there is a short range attractive potential. Magnons can 
be confined within the potential well (bound states) as 
long as the relative kinetic energy between the magnons 
is smaller than the interaction energy. 

This situation applies to the bound states in the AHC, 
the S = mode at all wavevectors and the 5 = 1 mode 
over limited wavevectors around the node points. These 
are positions where the kinetic energy is small compared 
to the binding potential. The 5* = 1 bound state is visible 
to neutron scattering and has a dispersion ||l^ 

iOBS = 2J- ^ (4cos2(fc(i/2) + 1) . 

The bound state is characterized by a probability am- 
plitude for the separation between the two magnons that 



TABLE IL Crystallographic data for Cu(N03)2 • 2.5H2O from 
reference jl5|. The room temperature lattice parameters are 
a = 16.453, b = 4.936, c = 15.963 A and /3 = 93.765° The 
Cu ions are at the 8f positions at (x,y,z) where x=0. 12613, 
y=0. 01352, z=0. 11376. The equivalent positions in the unit 
cell are: 



Atom x-pos y-pos z-pos 

1 X y z 

2 1-x l~y 1 - z 

3 1 — X y h ~ ^ 

4 X 1 — y \ + ^ 

5 }, + x i-y z 

6 2 ~ ^ 1 — z 

7 ^-x i-y i-2 

8 l; + x l + y \ + z 



drops exponentially with distance (see appendix for more 
details). These S = 1 bound states exist only over the 
range |n7r — kd\ ^ 7r/3, where n is an odd integer, see 
Fig. 1. 

MAGNETIC PROPERTIES OF CU(N03)2 2.5D2O 

The structural and magnetic properties of CN have 
been thoroughly investigated and shown to be near ideal. 
The structure of Cu(N03)2 • 2.5H2O was investigated 
by Garaj |l^ and Morosin and was shown to have 
a monoclinic crystal structure with space group /12/cl 
|l6| and the lattice parameters and crystallographic cop- 
per positions given in Table II. The deuterated form 
Cu(N03)2 -2.5020 we studied has low temperature (T = 
3 K) lattice parameters a— 16.1, h= 4.9, c = 15.8 A and 
/3 = 92.9° [0. 

The magnetism of CN arises from the Cu^"'' ions. Crys- 
tal electric fields from the oxygen ligands surrounding 
the Cu^"'' ions quench their orbital moments, leaving a 
near-isotropic spin-1/2 moment with g- values that show 
a small easy-axis anisotropy along the crystallographic b- 
direction; the values are g\\Y, = 2.33 and gj^b — 2.09 [p7| . 
Magnetic superexchange in this material is mediated by 
long double Cu-O-O-Cu exchange paths, which accounts 
for the rather weak exchange interaction observed in CN. 

The exchange couplings of CN were studied by Eckert 
et al. in |l^ and are illustrated in Fig. 3. The dominant 
magnetic exchange integral J is between pairs of spins 
(copper positions 7 & 8, and equivalent pairs, c.f. Fig. 3) 
forming dimers with a separation of 5.3 A. These dimers 
separated by crystal vector u = [ui, U2, M3] = uia-t-U2b-l- 
U3C are coupled together by exchanges J„; the only ex- 
change paths of appreciable strength are J,'i ,11,, con- 

nected via bonds between 1 & 7 and equivalent (bond 
length 6.2 A). This results in two sets of 5 = 1/2 al- 
ternating Heisenberg chains running in the [1.1,1] and 
[1,1,1] directions of the crystal with repetition every 
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d = [l,l,l]/2 and d' = [l,l,l]/2 respectively (repeat 
distance d = 11.3 A); the corresponding intradimer vec- 
tors are p = [0.252, ±0.027, 0.228] (where the x,y,z po- 
sitions from Table II have been used). Inelastic neutron 
scattering measurements |10 have recently further con- 
firmed the model alternating Heisenberg chain properties 
of CN and have shown that the dominant collective ex- 
citations are indeed the gapped triplet of magnons ex- 
pected for the AHC. 




FIG. 3: Positions of copper ions projected onto the ac plane 
for Cu(N03)2 -2.5020. The atom positions are those detailed 
in Table II. Two sets of identical chains run in the [1, 1, 1] and 
[l, 1, ll directions respectively. 



Approximate ground state energy 

Bulk magnetic measurements give information on the 
gap, exchange, and ground state energies of the spin 
chains in copper nitrate. Measurements of the effect of 
applied magnetic field show that spin flop (SF) ordering 
is induced in CN above a critical field B^x « 2.7 T, with a 
transition to full alignment at i?c2 ~ 4.3 T and there are 
no demagnetization effects in the zero temperature limit 
pof . As the orbital moment on the Cu^"*" ions is quenched 
by the crystal electric field and demagnetization effects 
are negligible, the field B^x can be used to directly give 
the excitation gap energy to the one-magnon states, and 
Diederix et al. report a value of A = 0.378 ± .007 meV 

The gap can be inferred from the transition field be- 
cause the magnons carry spin quantum numbers — 
1, 0, — 1 and are split into three dispersive modes shifted 
by a Zeeman energy with respect to each other. The 
ground state by virtue of its spin-0 quantum number is 
unaffected by the field and the transition occurs when 
the Zeeman energy of the lowest mode closes the gap 



and magnons condense into the ground state. The long 
range order itself is due to weak couplings between the 
chains. 

The high field transition, Bc2, yields further important 
information. It is where all the low-lying magnons are 
completely condensed into the ground state and the spins 
are fully aligned along the field. The fully aligned state 
is an exact eigenstate, and for an unfrustrated quantum 
magnet the transition field gives the sum of exchange 
couphngs 5/iB-Bc2 = J + T,u'^i = 0-580 ± .007 meV in 
the system |Q . An estimate of the ground state energy 
of CN can be made using these numbers. 

Using the low temperature isothermal magnetization 
M{B) = g^g {^^)b ^ ^ measure of the work re- 
quired to saturate the spin chains from the zero-field 
quantum ground state, an energy-per-spin can be in- 
ferred. The zero-field ground state energy-per-spin eo 
can be estimated via the formula eo ~ e/ — gn^SBc^ + 
M {B)dB where the fully-aligned energy-per-spm is 
e/ = S'V2-(J+Eu -^u) = 5Ms^c2/8 = 0.0725±.001 mcV 
for a 5* = 1/2 unfrustrated system. Utilizing the 270 mK 
data of Diederix et al. in Fig. 3 of ||l9| (measured using 
proton resonance) to determine the integral over magne- 
tization gives an experimental ground state energy-per- 
spin eo = -0.174±.004 meV. This is essentially the T = 
result, as the gap activation energy corresponds to 4.4 K. 

To estimate thermodynamic properties we approxi- 
mate the sum of the mterdimer exchanges by the sin- 
gle coupling aJ ~ V „ J[, of equation (1). Using the 
0{a^) expansions |l^ for A(a) and eo gives J = 0.455 ± 
.002 meV and a = 0.277 ± .006; in agreement with the 
resuhs of |l| and |§, J = 0.45 meV and a = 0.27. 
Our calculated values of the thermodynamic parame- 
ters J + Eu "^u = 0-581 meV, A = 0.379 meV and 
eo = —0.172 meV agree within error with the experi- 
mental values. 



EXPERIMENTAL METHOD 



Neutron scattering 



The inelastic neutron-scattering cross-section |£1| 

(5a/3-QaQ/3)5"''(Q,C^), 

is proportional to the dynamical response iS"^(Q,cl'), 
where Q is the wavevector transfer, F {Q) is the mag- 
netic form factor, N is the number of scattering centers, 
the constant amag = 0.2896 b, ki and kf are the mo- 
menta of initial and final neutron states respectively, g 
is the Lande (7-factor, and a = x,y, z are Cartesian coor- 
dinates. The dynamical response is the space and time 
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Fourier transform of the spin-spin correlation function 

where i and j labels sites of the system. For the AHC, 
equation (1), spin conservation and isotropy in spin 
space ensure that S"^{Q,lu) = for a 7^ /3, and all 
diagonal spin components are equivalent S^^{Q,uj) — 
Syy{Cl,u) = 5^^(Q,w). At T = this is given by 

S-^Q,u;) = ls+-iq,Lo) = |(vI.,(fc)|5+|*G)P^(c^ - 

A 

where A label the eigenstates of H and 




is the Fourier transformed spin creation operator. 

The action of the neutron is to flip a spin and so cre- 
ate a localized spin-1 polarization in the chain, and the 
strength of scattering to particular states is determined 
by their overlap with this spin flip state. This means 
that the multiparticle states will be sampled with parti- 
cles created close together. This is where interactions of 
the particle wavefunctions are most important, making 
neutrons a sensitive technique for looking at overlap ef- 
fects. An interesting consequence of this is that the short 
range interactions between particles can have a large in- 
fluence on measured correlation functions with little ef- 
fect on the thermodynamics of the spin chain, an effect 
noted for spinous in uniform chains |2^. Similarly, the 
thermodynamic influence of bound states at low temper- 
atures vanishes as 1/A'' yet they have a finite scattering 
cross section as discussed below. 

Neutron scattering measurements 

We made our measurements of the inelastic scatter- 
ing cross-section of CN using the SPINS cold neutron 
triple-axis spectrometer at the NIST Center for Neutron 
Research. The same high quality sample of copper ni- 
trate used by Xu et al. was utilized. This 14. Ig sam- 
ple consists of four coaligned single crystals of CN, with 
deuterium substituting for 92% of hydrogen. The substi- 
tution of D for H was made as it reduces significantly the 
background from incoherent scattering of neutrons but 
does not change the magnetic properties of the material. 
The sample was mounted with (/i, 0, Z) as the scattering 
plane in a pumped '^He cryostat at a base temperature 
of 300 mK. This temperature is an order of magnitude 
smaller than the gap energy (~ 4.4 K) and the collective 



quantum ground state is almost entirely free of thermally 
produced magnons, with a population of these numbering 
less than one per two million dimer sites. 

The two-magnon scattering is expected to be weak as 
the neutron matrix element to it is of 0{a'^) from the 
ground state so the spectrometer was set up in an open 
configuration to gain maximum scattered signal, and the 
only collimator included in the setup was of 80' between 
monochromator and sample. A vertically focused py- 
rolytic graphite PG(002) array monochromated the inci- 
dent neutrons (energy Ei, wavevector k^) and a horizon- 
tally focused array composed of eleven independently ro- 
^ tatable PG(002) blades was employed to analyze the scat- 
'\;e'red neutrons (i?/, k/). A cooled Be filter was placed in 
the incident beam before the sample to remove higher- 
order contamination from the beam. The actual neutron 
energy transfer to the sample being hw = Ei — Ef and 
the wavevector transfer is Q = k; — k/. 

Measurements of scattering cross-section were made 
by fixing the final energy at Ef = 2.5 meV {kf = 1.10 
A-i) and scanning incident energy Ei at various fixed 
wavenumber transfers along the chain, fc = Q • d (i.e. the 
component of the scattered wavevector along the impor- 
tant chain direction). Although there are actually two 
types of chain in CN (with repeats b = [1, 1, l]/2 and 
b' = [l,l,l]/2) this is not important in our case - we 
study the {h, 0, 1) scattering plane where the chains give 
identical contributions (see Fig. 3). 

With an open scattering configuration instrumental 
resolution is an important consideration. The spec- 
trometer resolution represents the spread in coordinate 
space (Q,w) sampled by the instrument at each mea- 
sured point. The energy resolution of the spectrome- 
ter is of Gaussian profile with a full- width-half-maximum 
(FWHM) at Ef = 2.5 meV and hw = 0.8 meV of ~ 0.10 
meV. The Q resolution is dominated by the wide angular 
acceptance (14°) of the analyzer on the scattered side - it 
is highly elongated along a direction within the scattering 
plane that is perpendicular to the scattered wavevector 
- and approximating the measured angular dependence 
by a Gaussian profile gives a FWHM of ~ 0.2 A-i. 
This resolution width is very considerable, however by 
using calculated scan trajectories that maintain the final 
wavevector kf along the crystallographic (1,0,1) direc- 
tion, so as to integrate over the nondispersive directions 
between chains, good wavenumber resolution in k along 
the important quantum spin chain directions (estimated 
at of order 0.02 A^^) is maintained. 

RESULTS AND ANALYSIS 

Fig. 4 shows some of the results of scans in energy 
performed on CN: Panel (a) shows a scan at the an- 
tiferromagnetic zone-center, k = 2iT/d, taken at T = 
300 mK. This is the wavenumber along the chain where 



FIG. 4: (a) Low temperature scattering &t k = 2-k / d. The 
dashed line is a fitted background and the sohd line is a fit to 
the scattering described in the text, (b) Two magnon scatter- 
ing with background subtracted off. The solid line is a fit (see 
text). The solid bar indicates the instrumental resolution, 
(c) Two- magnon scattering for k = Stt/c! with nonmagnetic 
background subtracted off. 

the magnon energy is a minimum, c./. Fig. 1. Strong 
elastic scattering from incoherent nuclear processes is 
clearly seen as well as a peak at 0.4 meV as expected 
for the one magnon mode |^| , close to the dimer excita- 
tion energy J = 0.45 meV. A second much weaker peak 
appears at roughly double the dimer energy at about 
0.9 meV, which is where two-magnon scattering is ex- 
pected. On heating up the sample both the 0.4 meV 
and 0.9 meV peaks disappear identifying these as being 
magnetic in origin. 

Heating up the sample also identifies the non-magnetic 
background contribution (dashed line in the figure) which 
consists of the incoherent nuclear peak, modelled by a 
Gaussian centered at zero-energy, and a broad contribu- 
tion from thermal diffuse scattering from the analyzer 
which is well-characterized by a power-times-Lorentzian 
(broad, quasielastic) component decaying from zero- 
energy. As the background is large compared to the two 
magnon signal it was studied in depth at different tem- 
peratures and wavevectors. 

Panels (b) and (c) of Fig. 4 show this peak with the 
modelled nonmagnetic background subtracted at A: = 
2-K / d and k — Ztt / d, respectively. This feature is con- 
siderably weaker than the one-magnon scattering. The 
FWHM of the peak narrows from 0.18(1) meV at 27r/d, 
panel (b), to 0.11(1) mcV at the zone-boundary (37r/(i), 
panel (c); behavior indeed consistent with two-magnon 
scattering. 

Fig. 1 shows the energy corresponding to a weighted 
average of scattering as a thick grey line. It is notable 
that at /c = 37r/d to a good approximation the neutrons 



couple only to the bound mode, so that nearly all the 
scattering weight is in it, not the continuum. The calcu- 
lated neutron scattering intensity from the bound state 
is ~ 2% of the one-magnon intensity which agrees with 
the data in Fig. 4. 

The one- and two-magnon scattering at 300 mK was 
scanned from /c = 7r/d to Stt/c? in steps of it /Ad. The 
background subtracted data are plotted in the upper 
panel of Fig. 5. The lower panel shows the calculated 
magnetic scattering based on the ID perturbation theory 
of the Appendix with the estimated parameters J = 0.45 
and a — 0.27. The calculation includes the correct dimer 
structure factor effects and uses the true scan trajecto- 
ries in conjunction with the dimer coordinates given in 
Table II. Corrections for the Cu^+ magnetic form factor 
|p3| have also been made. Instrumental line broadening 
has been included as well by convolving the theoretical 
calculations with the estimated instrumental resolution 
and sample mosaic spread. The calculation is directly 
comparable with the data in the upper panel of Fig. 3, 
and at a qualitative level there is good agreement with 
experiment. 

Next we consider the measured one- and two-magnon 
scattering in more detail and relate this to the physical 
picture presented by perturbation theory. A quantitative 
comparison between theory and data is shown in Fig. 6. 
The measured positions of one- and two-magnon peaks 
are plotted in the left panel. Energies, widths and in- 
tensities for each peak were extracted by least-squares 
fitting of Gaussians. In fact the wavevector and energy 
resolution was not sufficiently good in this experiment to 
distinguish details of line shape and the energies, widths 
and intensities represent the meaningful content of the 
measured data. We examine the one-magnon scattering 
first. 

One-magnon scattering 

Dispersion: The measured dispersion of the one- 
magnon states is shown in the left hand panel of Fig. 
6. Considerable dispersion of the one-magnon modes 
around the dimer energy (0.45 meV) is evident as ex- 
pected. Although the one magnon dispersion has been 
calculated to high order 0[a^) ||ll| previously, the small 
value of a in CN means that the one-magnon dispersion 
should be well approximated by the 0{a) result, equa- 
tion (2) . In fact the dispersions in CN measured by Xu et 
al. [0 and Stone et al. show that it is well described 
by 

u 

with J = 0.442(2) meV the dimer coupling, J[iii]/2 = 
0.106(2) meV the along chain coupling, plus additional 
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FIG. 5: (Color) Upper panel shows a color filled contour 
plot of the measured data with nonmagnetic background sub- 
tracted. Intensity is on a linear scale indicated by color, 
going from dark red (minimum) to light yellow (maximum). 
The two-magnon scattering has been enhanced by a factor of 
100 to make it visible on the same scale. Lower panel shows 
the calculated scattering using perturbation theory with cor- 
rections for instrumental resolution, multiple scattering and 
magnetic form factor. 

weak inter dimer couplings o o] — 0.012(2) meV and 
Jj'p p -^^2] — 0.018(2) meV. The alternation ratio that Xu 
et al. consider a = Jywm^l J = 0.240(5) is smaller than 
that found from the magnetization data discussed above 
which is presumably due to the neglect of interchain cou- 
pling effects in the analysis of the latter. 

The solid line through the one magnon dispersion in 
Fig. 6 is that calculated using the results of Q and it 
gives a reasonable account of the data. The small dis- 
crepancies from the peak-centers measured here are at- 
tributable due to the effects of the instrumental resolu- 
tion which averages over a large swathe of the interchain 
dispersion modulation. 

Intensity: The lower right-hand panel of Fig. 6 
shows the one-magnon intensities extracted from fit- 
ting. The neutron scattering matrix element (Q) = 



\{^x{k)\S^\^G)\'^ for excitation of the one-magnon states 
in the AHC to 0{a^) is 

St.n.giQ) = (1 - C0S(Q • p))- 

r (1 _ - ^a^) + (ia - - ^a^) cos(fcd)+ 1 
\ (^a^ + jga^) cos{2kd) + cos{3kd) j 

(4) 

The leading order scattering process is from the bare 
dimer component of the ground state, and the a/2 ■ 
cos{kd) component in the one magnon structure factor 
arises from an 0{a) two-dimer excitation in the ground 
state. The dynamical structure factor then reflects both 
the composition of the ground and excited states leading 
to a complex wavevector dependence. 

An interesting aspect of the one-magnon intensity 
5^7ag(Q) noted in |l^ is that the spin structure factor 
comes in only as a (1 — cos(Q • p)) modulation. So where 
Q • p = 27rn (n integer) the magnetic intensity should 
be zero. Although this situation occurs for the measure- 
ments at the wavenumber k = 3.97r/(i, the scattering 
does not go to zero because of residual intensity from 
secondary elastic scattering from incoherent processes. 
Such residual scattering was also observed in CN by Xu 
et al. [|l^ . The solid line in the figure then is the scatter- 
ing intensity predicted using equation (4) with J — 0.45 
meV and a — 0.27 including secondary scattering, as well 
as instrumental resolution and magnetic form factor. It 
is seen to account very well for the observed scattering 
intensity. We now consider the two-magnon scattering. 

Two-magnon scattering 

Integrated peak intensities, widths and positions ex- 
tracted using the fitting method described above for the 
two-magnon scattering are also plotted in Fig. 6. The 
peak positions are plotted in the left panel of Fig. 6 (filled 
circles). They are nearly dispersionless at an energy of 

0.86 meV. The peak widths are plotted in the top right- 
hand panel. These are largest around k = 2Tr/d and in/d 
where the continuum is expected to be widest, and are 
near resolution limited around 7r/d, STr/d, and bir/d. Also 
shown is the integrated peak intensity J dujl{uj), where 
is the intensity, which peaks around k = 27r/d. 

Intensity: The integrated peak intensity of the mea- 
sured two magnon scattering includes a sum over the 
two-particle continuum and = 1 bound mode. Tak- 
ing account of the density of momentum states with en- 
ergy transfer, the scattering intensities for the continuum 
(eqn. (17)) and bound (eqn. (18)) states as described 
in the Appendix were computed. As interchain effects 
are effectively integrated over in the two-magnon scat- 
tering, and this results in line broadening rather than 
shifts in energy; we thus have implicitly included inter- 
chain coupling effects in our definition of a — 0.27 for (1). 
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The energy integrated two-magnon intensity, / du!l{uj), 
(right middle panel) shows a more complicated k depen- 
dence than the one magnon scattering. The compari- 
son with the 0{a'^) calculation looks qualitatively simi- 
lar to the data, however it underestimates the scattering 
at fc « 97r/2d and overestimates it at 2Tr/d, which may 
indicate that higher order terms in the scattering ampli- 
tude are important. It is notable that the two-magnon 
intensity is very strongly dependent on the spatial ar- 
rangement of magnetic ions. 

Center: the fitted peak centers are compared with the 
computed weighted center (uj) = f dujl{uj) x uj (gray 
band) for the 0{a^) perturbation theory in the left panel 
of Fig. 6. The nearly dispersionless extracted positions 
(grey filled circles) are located at the calculated weighted- 
average energies (grey band) replotted from Fig. 1. The 
fact that the weighted center lies below the center of the 
continuum is a direct result of the movement of scattering 
weight towards the lower boundary due to the magnon- 
magnon interaction. 

Width: the peak widths obtained from the fits are 
shown in the top right-hand panel. The solid line rep- 
resents that calculated using the perturbation theory. It 
is the sum-in-quadrature of the instrumental resolution 
width in energy and the variance a of the theoretical in- 
tensity where tr^ = / d(jjl{uj)x{uj—{uj))'^ . The calculation 
is seen to provide a good account of the data. 

Bound mode: One of the most interesting aspects of 
the multiparticle states in the AHC is the existence of the 
bound mode below the two-magnon continuum. The pre- 
dicted wavevector dependence of the intensity (see Ap- 
pendix) is 

'5^s"(Q) = (f)'[l-4cos2(fed/2)] X 

[sin (Q • (p + d)/2) + 3 sin (Q • (p - d)/2)]' . 



and should be visible around k — Sn/d. The binding en- 
ergy of the S* = 1 state, neglecting interchain coupling, 
is predicted to he ^ Eg ^ J il^ ~ = 0-01'' 

meV for CN. The scattering in CN around k = Sn/d 
is centered at 0.852 ± .007 meV, which gives a binding 
energy of Eb = 0.03 ± .02 meV. In addition the scat- 
tering is near resolution limited, as expected for a well- 
defined mode. However, although the energy and inten- 
sity around k — Sn/d lend support to binding around 
this bandwidth minimum, the experimental error means 
this does not constitute definitive proof of the eifect in 
CN. 



practical techniques for calculating the multiparticle exci- 
tation spectra and cross sections for realistic spin models. 
Our perturbation theory, although useful for interpreting 
results, is of too low order to quantitatively account for 
our measurements and in addition does not include inter- 
chain coupling. Very powerful linked-cluster-expansion 
techniques have recently been introduced that allow mul- 
tiparticle spectra to be calculated to high order |2^ 
for the AHC and the extension of these to calculations 
of the neutron scattering cross-section would be a signifi- 
cant development. Analytical approaches based on Green 
function techniques may also prove fruitful. Secondly, 
measurement of weak multiparticle signals preludes the 
access of neutron scattering to measuring bound two- 
magnon states. Interactions such as next nearest neigh- 
bor coupling should further stabilize bound modes and 
make measurement of these easier. A question which re- 
quires further investigation is the stability of such bound 
modes to thermal fiuctuations and also interchain cou- 
pling. An alternative route for investigating the phe- 
nomenology of particle binding is through solitonic sys- 
tems such as the ID Ising chain with small XY-like terms 
prf . In this case the coupling of neutrons to pairs of 
5 = 1/2 solitons (also called spinous) is at zeroth or- 
der and therefore strong. Binding of solitons only occurs 
when extra terms are included in the Hamiltonian, such 
as exchange mixing |2^, next-nearest neighbor coupling 
p9[, or transverse field |30 . Evidence for this binding 



DISCUSSION 

Our measurements establish the feasibility of study- 
ing weak multi-magnon states using neutron scattering 
and raise a number of important issues. Firstly the in- 
troduction of experimental data highlights the need for 



phenomenon in and Ising-like chain has been observed 
recently in neutron scattering experiments |Q . 

Time-of-flight (TOF) neutron spectrometers give po- 
tentially much better energy resolution than conventional 
triple-axis instrumentation and could provide definitive 
proof of the bound state in the alternating chain by re- 
solving the bound mode from the continuum. Previously, 
TOF techniques have proven successful in the study of 
similar binding effects at the bandwidth minimum of the 
two-soliton continuum scattering of the 5* = 1/2 XXZ 
Ising chain material CsCoCla |2^] . Although limited neu- 
tron fluxes may make such measurements difficult for CN 
|l0| these should be feasible and we plan to make such 
measurements in the near future. 



SUMMARY 



In summary, we have used inelastic neutron scattering 
to investigate the ground and excited states of the near- 
ideal alternating Heisenberg chain material Cu(N03)2 • 
2.5D2O, and also derived the scattering analytically to 
lowest order in perturbation theory. Our measurements 
are consistent with the predictions of this model for sev- 
eral magnetic properties of this system, including the 
ground state energy, one- and two-magnon excitation 
spectra and intensities, and possibly the existence of a 



As total spin St = E^li^™, and = E™=i ^™ 
are constants-of-the-motion for the Hamiltonian H they 
organize the Hilbert space. For notation we introduce 
dimer creation operators a'^\Gm) = \o-m) where a — 
1,0,1 label excited dimer states. In real space we de- 
note the singly excited states with quantum numbers 
{St,S^) = {I, a) as |a),„ — a^^^o and the doubly ex- 
cited with \{St, S^))m,i^ as the Clebsh-Gordan combi- 
nations of excitations at dimer sites m and m + v i.e. 
1(0,0))™,. = l/\/3{l+T++, - 0+0++, + T+ 1++JVI/0, 

1(1,1))™,, = 1/V2{0+1++, - l+0++J*o, etc. 

As the alternating chain has translational symmetry 
plane-wave states prove convenient 



\a)k 



1 



(6) 



FIG. 6: Comparison of theory and data. Left panel shows 
fitted positions of observed sattering using perturbations the- 
ory (see text). Right lower panel shows fitted one magnon in- 
tensity compared with perturbation theory (see text). Right 
upper panel shows a comparison of the two-magnon intensity 
with perturbation theory (see text). 



two-magnon bound state. Much experimental work re- 
mains to be done to establish the phenomenology of bind- 
ing in isotropic ID systems. 
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APPENDIX 



Ground state 



and 



\iS,Sn)k.. ^ 4f E e™^'l(^,^^))™,.- (7) 

where the allowed momenta are k„ = 2m:/Ndd, where 
n are integers from —Nd/2 to Nd/2. The action of the 
Hamiltonian on the basis states has been considered in 
|pd| and the ground state to 0{a) is 



%[vI/o-a^^ 1(0,0)), 

m — 1 



(8) 



where ^ 1 ~ {3/128)a^ Nd. 

Neutron scattering measures the square of 
the expectation value of the spin operator 
S+{Q) = (2A^d)-^/'E™=iEp=±exp(zQ • r™,p)5+p 
between eigenstates, where Q is the wavevector transfer 
of the neutron. The action of this spin operator applied 
to the ground state is 

S+{Q)\^G)^Aci\l)k + Bci\{l,l))k, 



A 



Q 



1 

V2 



-iQp/2 



,-tQp/2\ 



) (1 + f cos(fcd)) , (9) 
) (e''^''* - 1) . 



where fc = Q • d is the wavenumber of the states ex- 
cited by this operator. The neutron scattering matrix 
element at T = to state A, is given by (Q) = 



Here we apply the perturbation theory of to two- 
magnon scattering. The single dimer eigenstates (la- 
belled G, 1, 0, 1) are listed in Table I and the ground state 
of the uncoupled dimers {a = 0) is a direct product of 
dimer ground states 



^-0 = 



n 



\G„ 



En 



3JNd 



(5) 



One-magnon states 

The one-magnon wavefunctions to order 0{a) are pd| | 
*imag = \a)k + ^{e^''' + im,a))k,i~ (10) 



y\/3 



J2 S„^^rn'e'"""'al\{0,0))rn,,^ 



m,m' — 1 
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and form an S' = 1 triplet with an energy gap A = J 
aJ/2 above the ground state and dispersion 



f-^imagCfc) = J — aJ/2 cos{kd). 



(11) 



Apphcation of the S'+(Q) operator gives the intensity to 
Oia) 

5i+„rag(Q) = (1 - C0S(Q • p))- 

{l + ^acos{kd)] . 



Two-magnon states 

Ignoring states higher than two-excited-dimer (as they 
do not contribute to the neutron-scattering matrix ele- 
ment to lowest order in a), the action of the Hamiltonian 

is fini 



H\a)k = ak\a)k + V2jl\{l,a)) 



k,l 



and 



1 /3|(1, a))k., + 7fc|(l, a))k.,+i + a))fc,.-i, ^'^t' 



I V27fc|a)fc + (/3-e)|(l,a))fe.i+7fcl(l,a))M, ^ = [12) 



where = J — {aJ/2) ■ cos{kd), (3 = 2 J, e = a J/4, 
jf, = — (1 -I- e"*'^''), and f denotes complex conju- 
gation. The excitation spectrum can be calculated by 
direct diagonalization of a large number of dimers and 
application of the matrix element above or by analytical 
solution. 

Approximate analytical wavefunctions for the S — 1 
states can be calculated using elementary scattering the- 
ory, see Ignoring the coupling to the one-excited- 
dimer states for the time being, the two-magnon wave- 
functions are 



*2mag(fc) = &^exp(i6ll/)|(l,l))fe,,. 



(13) 



where for 9 = ~i'^og{\J l\hk) — fcd/2 the time- 
independent Schrodinger equation reduces to solving the 
real and symmetric system of equations 



\hi = {(3- e)6i + 7^62 



(14) 



\b 



= {[3- e)bNd-i + Ik^Nd- 



where 7^ = \Jlklk — a J/2 • |cos(fcc?/2)| and the term 
exp(i0i/) serves to transform to the center-of-momcmtum 
(center R) frame where magnons are at = i? — v /2 and 



R-\- V 12 with total momentum K — ki 



kd. 



Magnon-pair-state solutions 

Magnon-pair-states comprise particles that are free 
at large distances and for particle conservation in one- 
dimension the interactions introduce a phase factor (j) on 
scattering. This state corresponds to 

K = ^0 (exp {ip^,iy) - exp {-i {p^iy ~ (Z-^))) (15) 

where the normalization constant Xq ~ 1/V7V^, is the 
relative momentum, and fi = 1,2, ...,iVjj — 1 index the 
eigenstates. The phases and momenta are determined 
by the boundary conditions of the particles and their 
interaction energy. The standard method of solving for 
these boundary conditions ^ is to introduce the 
single site coefficients 60 and Bn and set ebi — 7^.60 and 
e6jv^_i = ^}ANa and substitute into equation (14). The 



momentum p^ and phase 
60 and fo^r when 



solve for the constraints on 



7fe + eexp (-^^vj^j 



Ik + eexp 



(•- 



,Pli 



Nd 



(16) 



and the eigenvalues are \^ — 2J — aJ cos{kd) cos{p^) + 
0{1/Nd). In the time-dependent Schrodinger picture this 
state corresponds to two particles with wavepackets at 
positions = R + v /2 and rj — R ~ v /2 and with mo- 
menta fci — {K +p)/2 and ^2 ~ [K — p)/2 that scatter 
via the S-matrix ^kiM — ~ 6xp(i0j.^_j.^), and the eigen- 
spectrum is 



^kiM — ^lmag(fcl) + 'i^lmag(fc2) 



(17) 



with a;imag(^) = J ~ a J/2 • cos(fcd) as above. 

Including the coupling to the |l)fe states as a pertur- 
bation gives the approximate wavefunction 



2mag 



{k,y) 



ac 



2V2 



^(1 + ' 



-ikd' 



\^)k+ ci^i(i'i))fe,- 

(18) 



where — b^ exp [ikdv /2). The neutron scattering ma- 
trix element can then be computed straightforwardly for 
a large Nd system by evaluating the closure with equation 
(8). 



Bound state solutions 

When 7j. < e, (j)^ cannot be solved with /i = 1 and Nd— 
1, and the c\, and c^'^~^ solutions comprise exponentially 
decaying bound states solutions below the two magnon 
continuum. The 5=1 bound state wavefunctions and 
energy have previously been given by Uhrig and Schulz 
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|g[ , and Damle and Nagler ||lJ 
(18) with 



It is specified by equation 



,1 h"^ ^ I'k , \ 

V 7fc 



(19) 



where exp(— k) = ^7fe/e which has dispersion energy 
UJBS = 2J- ^ (4cos^(fcd/2) + 1) . 

An interesting feature of the bound state solution is that 
it only exists over the range \mT — kd\ ^ 7r/3, where n 
is an odd integer, and so there is an S* = 1 bound mode 
only for small wavevectors around the narrowest part of 
the continuum. Using our wavefunction we calculate the 
neutron scattering strength from the bound state to be 

'5^s"(Q) = (f)'[l-4cos2(fcd/2)] X 

[sin (Q • (p + d)/2) + 3 sin (Q • (p - d)/2)f . 
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